Majorana solution of the Thomas-Fermi equation 
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We report on an original method, due to Majorana, leading to a semi-analytical series solution of the 
Thomas-Fermi equation, with appropriate boundary conditions, in terms of only one quadrature. 
We also deduce a general formula for such a solution which avoids numerical integration, but is 
expressed in terms of the roots of a given polynomial equation. 
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In 1928 Majorana found an interesting semi-analytical 
solution of the Thomas-Fermi equation (!]] which, unfor- 
tunately, remained unpublished and unknown until now 
(see |0 ) . We set forth here a concise study of such a so- 
lution in view of its potential relevance in atomic physics 
as well as in nuclear physics (as, for example, in some 
questions related to nuclear matter in neutron stars 0). 
The problem is to find the Thomas-Fermi function fix) 
obeying the differential equation: 
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with the boundary conditions: 

<f(P) = 1 
f(oo) = 



(1) 



(2) 
(3) 



An exact particular solution of Eq. (|l|) satisfying, how- 
ever, only the condition (||), was discovered by Sommer- 
fold 19: 
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(4) 



This can be regarded as an asymptotic expansion of the 
desired solution, and Sommcrfeld himself considered a 
"correction" (in some sense) to Eq. in such a way 
to take into account the condition (0)^ However, this 
"corrected" approximate solution had a divergent first 
derivative for x = Q. 

As will be clear below, Majorana solution can be consid- 
ered as a modification of Eq. ([!]) as well, but the method 
followed by him is extremely original and very different 
from the one used by Sommerfeld. 

Let us consider solutions of the Thomas-Fermi equation 
(pi) which are expressed in parametric form: 



x = x(t) 

If = <f(t) 



(5) 



To be definite, throughout this paper we will use a prime 
' or a dot ' to denote derivatives with respect to x or t, 
respectively. The strategy adopted by Majorana is to 
perform a double change of variables: 



ip(x) 



t , u(t) 



(6) 



where the novel unknown function is u(t). The relation 
connecting the two sets of variables (assumed to be in- 
vertible) has a differential nature, that is: 



t = t(x, if) 
u = u(ip, f') 



In such a way the second order differential equation (|i|) 
for if is transformed into a first order equation for u. 
Note, however, that in general Eqs. (0) are implicit equa- 
tions for t and u, since x and <f depend on them (one is 
looking for parametric solutions in terms of the parame- 
ter t and the unknown function u). For the specific case 
of the Thomas-Fermi equation, Majorana introduced the 
following transformation: 



t = U4- 1 ^x 1 / 2 ip 1 / 6 



(8) 
(9) 



Observe that Eq. (g) is reminiscent of the Sommerfeld 
solution, since it can be cast into the form: 
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^ x 3 



t° 



(10) 



The differential equation for u(t) is obtained by taking 
the t— derivative of Eq. (ph: 



du 
~dt 
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3 ip 
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(11) 



and inserting Eq. ([!]): 



1 



du 



1G 



1/3 



X ip 



-4/3 



4^ 

3 ip 



•r 



.3/2 



,1/2 



(12) 



By using Eq. (0) and Eq. (||) to eliminate X 1 / 2 and ip' 2 , 
respectively, we obtain: 



du 
~db 



4\ 1/3 tu 2 



■ x ip 



1/3 



(13) 



We have now to express the quantity xtp 1 ^ in terms of 
t,u. From Eq. (§), 



2 = 144 1 / 3 £ 2 V3 1 / 3 , 
by taking the explicit t— derivative of both sides, 



= 144 1/3 



2tip~ 1/3 + t 2 



after some algebra we get: 

xtp 1 / 3 = 144 1/3 



2t 



l-t 2 i 



(14) 



(15) 



(16) 



By inserting this result into Eq. (|13|) we finally have the 
differential equation for u(t): 



du 
~db 



tu 2 - 1 
1 - t 2 u 



The condition (||) implies, from Eqs. (| 
for x — and: 



«(0) 
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(17) 
that t = 

(18) 



where <p' = ip'(x = 0). The initial condition to be satis- 
fied by u(t) for the univocal solution of Eq. ( |l7j ) is ob- 
tained from the boundary condition (0) by inserting the 
Sommerfeld asymptotic expansion (|J) into Eqs. (||),(|j]). 
For x — > oo we have t — 1 and 



u(l) 



1 



(19) 



We then easily recognize that the branch of u(t) giving 
the Thomas-Fermi function (in parametric form) is the 
one between t = and t = 1. In this interval we look for 
the solution of Eq. ( |l7| ) by using a series expansion in 
powers of the variable r = 1 — t: 



a Q + a\T + a 2 r 2 + a 3 t 3 + . . 



(20) 



From the condition ( |l9| ) we immediately have: 

oo = 1 . (21) 

The other coefficients are obtained by an iterative for- 
mula coming from the substitution of (BOh into Eq. (17): 



k+l _ 



(22) 



where: 

A(k, I) = a k [(I + l)a l+1 - 2(1 + 4)a ; + (I + 7)a;_i] + 

- (k + I + l)Si a k+ i +1 + SS k0 6i (23) 



(we define a„i = 0). Eq. 
form (k + I = m, I = n): 



can also be cast in the 




A(m — n,n) J r 1. 



(24) 



so that, for fixed m, the relation determining the series 
coefficients is the following: 



E 

n=0 



_„ [(n + l)a„+i - 2(n + 4)a„+ 



+ (n + 7)(1 - 8 n0 )a n 



(m + l)a m+ i 



JmO 



(25) 



(we have explicitly used that a_i = 0), with m 
0, 1, 2,3,.. .. The equation ((H) for m = 0: 



(oo - 1) [ox - 8 (oo + 1)] = 

is identically satisfied due to Eq. (|2l|). For m 
have a second degree algebraic equation for a±: 



18 ai + 8 







(26) 
1 we 

(27) 



of which we have to choose the smallest root (we are 
performing a perturbative expansion): 



Ol 



73 



(28) 



The remaining coefficients arc determined, using Eqs. 
( pl| ) and fl28|) , by linear relations. In fact excluding the 
cases with m = 0, 1, after some algebra Eq. ( |25| ) can be 
written as: 



1 



'm-2 

E 



a m -« [(n + l)t 



n+l 



+ 



2(m + 8) - (m + l)ai 
— 2(n + 4)o„ + (n + 7)o„_i] + a m _i [(m + 7)+ 
-2(m + 3)ai]+a m _ 2 [(m + 6)oi]} . (29) 

Note that the sum in the RHS involves coefficient aj with 



indices i < to— 1, so that the relation in (29) gives explic- 
itly the value of a m once the previous to — 1 coefficients 
a m _i, a TO _2, ...,a2)Oi (and oq) are known. 
The series expansion in ( p0[ ) is uniformly convergent in 
the interval [0, 1] for r, since the series made of the coef- 
ficients only, X^^Lo a ™' ^ s convergent. In fact, by setting 
r = 1 (t = 0) into @, from Eq. (ph we have: 



^ = I 16 



1/3 oo 

E 

n=0 



(30) 



fc=0 ;=o 



which shows that the sum of such a series is deter- 
mined by the (finite) value of ip' (ip' Q ~ 1.588 and thus 
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TABLE I. Numerical values for the first 20 coefficients for 
the series expansion of the function u(t) in Eq. (0). 



ai 


0.455996 


an 


0.0316498 


n-2 


0.304455 


012 


0.0252839 


az 


0.222180 




0.0202322 


04 


0.168213 


ai4 


0.0162136 


O5 


0.129804 




0.0130101 


«<; 


0.101300 


a iii 


0.0104518 


07 


0.0796352 


an 


0.00840559 


as 


0.0629230 


ais 


0.00676661 


a<> 


0.0499053 


aig 


0.00545216 


«io 


0.0396962 


a2o 


0.00439678 



S^Lo a ™ — 2.7746). Note also that the coefficients a n 
are positive definite and that the series in (|3^) exhibits 
geometric convergence with a n /a„_i ~ 4/5 for n — > 00. 
The numeric values of the first 20 coefficients are reported 
in Table I. Given the function u(t) we have now to look 
for the parametric solution x — x(t), tp — tp(t) of the 
Thomas-Fermi equation. To this end let us put: 



ip(t) = cxp 



w{t) dt 



(31) 



where w(t) is an auxiliary function to be determined in 
terms of u(t), and the condition (^) (or ip(t = 0) = 1) is 
automatically satisfied. By inserting Eq. (|3l]) into Eq. 
(^) and using ( |l6| ) we immediately find: 



6ut 



l-t 2 u 



(32) 



Summing up, the parametric solution of Eq. ([!]) with the 
boundary conditions (|^), (|^) takes the form: 



with: 



x(t) 

At) 




tit 



1 - t 2 ; 



■ dt 



(33) 



(34) 



and u(t) is given by the series expansion in (20) with 
the coefficients determined by (|2l|), ( pq ) and (29|). Eq. 
( |33| ) represents the celebrated Majorana solution of the 
Thomas-Fcrmi equation; it is given in terms of only one 
quadrature 0. 

We have performed numerically the integration in ( |34| ) 
stopping the series expansion in (pOj) at the terms with 



n = 10 and n = 20, respectively, and compared the para- 
metric solutions thus obtained from (]3^) with the exact 
(numerical) solution of the Thomas-Fermi equation. We 
found that the two Majorana solutions approximate (for 
excess) the exact solution with relative errors of the order 
of 0.1% and 0.01%, respectively. 

We can also obtain an approximate (by defect) analytic 
solution by inserting the series expansion (^p| ) into the 
expression (F}4]): 



I(t) = 



u(l - t) 
_ t l-(l-r) 2 u 

1 + biT + b 2 r 2 



dr = 



CiT + C 2 T Z + 



dr 



with: 



bn &TI 



('n b n —i b n 
for n > 1, while bo = 1 and cq = 0. Note that: 
b„ < for n > 1 



< 



for n > 1 



(35) 



(36) 



(37) 



(and bo, a > 0). If we neglect 0(t ) terms in (p5[), the 
quantity lit) is approximated by: 

l(t) 



-t- - log(l-t) 



(38) 



Ci Ci 

and, in terms of the original a„ coefficients, the approxi- 
mate parametric solution of the Thomas-Fermi equation 
is: 



x(t) 



Suit 2 (i-ty 



(39) 

tp(t) = (l-t)^e 6 ^T' . 

In Fig. we compare the above solution with the exact 
(numerical) one. 

More in general, we can truncate the series in ( p5| ) to 
a certain power r fe and thus the integrand function is 
approximated by a rational function: 



1 + bir + b 2 r 2 



C\T + C 2 T* 



b k r k _ P(r) 
" Q(r) 



CkTk 



(40) 



Let us then assume that the roots t% (i = 1, 2, . . . , k) of 
the polynomial in the denominator, 



Q(n) 







(41) 



are known, so that we can decompose the function F[t) 
in a sum of simple rational functions ^ 



*Eq. (|3^) is, probably, the major result of this paper and was 
obtained by Majorana. What follows is, instead, an original 
further elaboration of the material presented above 



^For simplicity we are also assuming that all the zeros of 
Q(t) are simple roots, as it is likely in the present case. How- 
ever, the generalization to the case in which multiple roots 
are present is straightforward. 
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FIG. 1. The Thomas-Fermi function <p(x). The thin (up- 
per) line refers to the exact (numerical) solution of Eq. (|l|), 
while the thick (lower) one corresponds to the parametric so- 
lution reported in Eqs. (p9|). 



1 



Q(r) 
and: 



Ck 



n 



T-T 2 



T-Tk 



(42) 



F(t) = — (l + bir + b 2 r 2 

Ck 



T-Tl 



T -T2 



b k r k ) 

fk 



The expressions for the coefficients /j (i = 1, 2, . 
terms of the roots t* are as follows: 



k 

n 

i = i 

l ^ i 



n 



By inserting the decomposition ( (43[ ) into Eq. ( |35| ) , 
tegral T(t) is thus given by a (double) sum whose 
element has the following form: 



T n 1 — Ti — t 

dr = -rf log 



(43) 
, k) in 

(44) 



the in- 
generic 



l-n 



+ £(? j^dl-rj-il-n-t) 1 ) . (45) 

Then, in general, the parametric solution of Eq. ([j]) can 
be formally written as: 



n—l 



x{t) = Vult 2 e 2 {t) 
<p{t) - e-\t) 
where e(t) — e x ^ is approximated by: 



(46) 



e(t) 



n ^ 

i—l 

n 

nnn 

i=l n =l 1=1 



1 — Ti — t\ c " 



(47) 



Obviously, using the method described above, the exact 
result is recovered in the limit k — * oo. This procedure 
can, however, be employed for getting approximate but 
accurate solutions of the Thomas-Fermi equation since, 
as it is clear from above, we have translated a numer- 
ical integration problem (see Eq. ([34])) into the one of 
a numerical search for the roots of the polynomial Q(t). 
Note also that we already know one of such roots (namely, 
Ti = 0) given the particular form of Q(t). This implies 
that, since the general solution of a fourth-degree poly- 
nomial equation in terms of radicals is known, from ( |4^ ) 
and (E7l) we can get an analytic approximate solution by 
considering terms in the series in (35) up to order 0(t 5 ), 
thus obtaining a certainly much better approximation to 
Eq. @ than Eq. @. We do not report here the 
explicit form of such a solution because of its very long 
expression. 

Summarizing, in this paper we have reported on an 
original method, due to Majorana, forwarding a semi- 
analytical solution of the Thomas-Fermi equation (|l|) 
with boundary conditions (^), (|J). The procedure ap- 
plies as well to different boundary conditions, although 
the constraint ^) is always automatically satisfied. This 
corresponds to physical situations present in atomic as 
well as in nuclear physics. We have further studied the 
Majorana series solution thus obtaining a general formula 
whose degree of approximation is limited by the one for 
searching roots of a given polynomial rather than to the 
one for integrating a rational function. 
The method used by Majorana for solving the Thomas- 
Fermi equation can be generalized in order to study a 
large class of ordinary differential equations, but this will 
be discussed elsewhere. 
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